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Abstract 



We study thermodynamic properties of the multiquark nuclear matter. The de- 
pendence of the equation of state on the colour charges is explored both analytically 
and numerically in the limits where the baryon density is small and large at fixed 
temperature between the gluon deconfinement and chiral symmetry restoration. The 
gravitational stability of the hypothetical multiquark stars are discussed using the 
Tolman-Oppenheimer-Volkoff equation. Since the equations of state of the multiquarks 
can be well approximated by different power laws for small and large density, the con- 
tent of the multiquark stars has the core and crust structure. We found that most 
of the mass of the star comes from the crust region where the density is relatively 
small. The mass limit of the multiquark star is determined as well as its relation to 
the star radius. For typical energy density scale of 10 GeV/fm 3 , the converging mass 
and radius of the hypothetical multiquark star in the limit of large central density are 
approximately 2.6 — 3.9 solar mass and 15-27 km. The adiabatic index and sound speed 
distributions of the multiquark matter in the star are also calculated and discussed. 
The sound speed never exceeds the speed of light and the multiquark matters are thus 
compressible even at high density and pressure. 
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1 Introduction 



All of the high energy experiments which fail to produce a free quark are strong evidences 
that the coupling constant of the strong interaction becomes nonperturbatively large at low 
energy and large distance. Quarks and gluons are said to be confined within hadrons and the 
colourless condition becomes a requirement of an assembly of quarks at low energy. However, 
when the energy or temperature scale of a system of quarks and gluons increases, the coupling 
of the strong interaction tends to be weaker and finally we expect the deconfinement to occur. 
In addition, if the quarks and gluons are compressed extremely tightly together, quarks could 
interact with neighbouring quarks and gluons equally and become effectively deconfined from 
the mesonic or baryonic bound state. In the latter case, the coupling could still be strong 
despite of the deconfinement. Nevertheless, we could also have the situation where gluons 
are deconfined but the quarks are not completely free due to the remaining Coulomb-type 
potential from gluon exchanges between quarks. 

Recently, the experimental results from collision of heavy ions suggested that the nuclear 
deconfinement phase might have been created in the laboratory and we might have pro- 
duced the quark-gluon plasma (QGP). The RHIC experiment revealed that the produced 
QGP behaves like fluid with very small viscosity. However, this property of small viscosity 
fluid is hard to be understood in the picture of QGP as the gas of free quarks and gluons. 
Additionally, lattice simulations show that QGP has relatively high pressure right above the 
deconfinement temperature T c which is again difficult to explain using the weakly coupled 
quarks and gluons gas [U El E]- It is possible that various coloured and colour- singlet bound 
states of quarks and gluons could exist in the plasma at the temperature (1 — 3)T C [U EJ [5]. 
The existence of the coloured bound states could explain the problems of high pressure, 
small viscosity, and the jet quenching of the QGP at once. 

Due to the large coupling of the strong interaction at low energies, a perturbative method 
has limited applicability to the high energy processes and phenomena. The development of 
the holographic principle and AdS / CFT correspondence [6] provides us with a new method to 
investigate the physics of strongly coupled nuclear matter both in the low energy regime and 
in the energy scale close to the deconfinement temperature. Holographic models of meson 
were proposed by Juan Maldacena, Soo-Jong Rey, Stefan Theisen, Jung-Tay Yee [3 EJ [9]. 
The Coulomb potential plus screening effect of quark and antiquark are calculated from the 
Nambu-Goto action of the string in the bulk spacetime at zero and finite temperature. For 
baryons, Witten, Gross and Ooguri [101 [TTJ proposed a holographic baryon to be a D-brane 
wrapping internal subspace of the background spacetime with N c strings connected and 
stretching out to the boundary. For AdS 5 x S 5 , the baryon vertex is a D5-brane wrapping 
the S 5 . The basic requirement is that a total of N c charges from the endpoint of the strings 
cancel with the charge of the vertex itself. A generalization of this condition allows more 
strings to go in and come out of the vertex, as long as the total charges from all of the string 
endpoints add up to N c [TJl [131 EH1 HI E] • Baryon vertex plus strings configuration in 
this case represent the holographic multiquark states. Generically they have colour charges 
but because of the confinement, they can only exist in the deconfined phase. 

The coloured multiquark phase can be studied in the general Sakai-Sugimoto model (SS) [TSl 
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\W\ in the intermediate temperature above the gluon deconfinement but below the chiral 
symmetry restoration temperature [2D]. It was found that the multiquark phase is ther- 
modynamically stable and preferred over the other phases in the gluon-deconfined plasma 
provided that the density is sufficiently large [17]. The situation of high density and moder- 
ate temperature could exist inside certain classes of compact stars and it is thus interesting 
to investigate the thermodynamical properties of the multiquark nuclear matter as well as 
their contributions to the stability of the dense stars. In this article, we will consider the 
hypothetical multiquark star which obeys the equation of state derived from the holographic 
multiquarks in the SS model. With the power-law approximation of the equations of state, we 
study its gravitational stability using the Tolman-Oppenheimer-Volkoff equation (TOV)[2Tj. 
The mass, density and pressure distributions are obtained numerically. The mass-radius 
relation and the mass limit are also discussed. Corresponding hydrodynamical properties 
such as the sound speed of the multiquark nuclear matter are explored within the star. The 
multiquark matters are found to be compressible throughout the entire multiquark star. 

This article is organized as the following. Section 2 describes the holographic setup for 
the multiquarks and the multiquark phase in the gluon-deconfined SS model. The thermo- 
dynamic relations and the equations of state of the multiquark nuclear matter are calculated 
and discussed in Section 3 and 4. In Section 5, the Einstein field equation for the spheri- 
cally symmetric star is solved to obtain the TOV equation. Assuming the equations of state 
derived in Section 3 and 4 for the multiquark nuclear matter, we explore the gravitational 
physics of a hypothetical multiquark star. A mass-radius relation is derived and some dis- 
cussion on the more realistic situation is commented. The adiabatic index and the sound 
speed of the multiquark nuclear matter within the star are studied. Section 6 concludes the 
article. 

2 Holographic multiquark configuration 

Since string theories in the bulk spacetime correspond to certain gauge theories on the 
boundary of that space, it is natural to find construction of the bound states of quarks in 
the form of strings and branes. While the meson is proposed to be the string hanging in the 
bulk with both ends locating at the boundary of the AdS space [7J, the baryon is proposed to 
be the Dp-brane wrapped on the S p with N c strings attached and extending to the boundary 
of the bulk space PHEE]. 

On the gauge theory side, hadrons exist in the confined phase as a result of the linear 
part of the binding potential. However, the bound states of quarks can actually exist in the 
deconfined phase at the intermediate temperatures above the deconfinement as well. Even 
though gluons are free to propagate and the linear potential is absent, the quarks can form 
bound state through the remaining Coulomb-type potential due to the colour charges of the 
quarks. 

The holographic model of non-singlet bound state was also proposed. As is demonstrated 
in Ref. [IT], we can modify the Witten's baryon vertex by attaching more strings to the vertex 
provided that the total number of charges of all of the strings are preserved to N c . Some 
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strings may extend along radial direction of the AdS space down to the horizon and some 
can extend to the boundary. We define the number of strings that extend to the boundary to 
be kh and the number of strings extending radially to the horizon to be k r . The restriction of 
kh and k r is due to the force condition of the string configuration (see Ref. [17] for details). 

In this article, we consider the holographic model of multiquarks in the Sakai-Sugimoto 
(SS) model [TSj [T§] similar to the configurations considered in Ref. [T7]. The background 
metric of the bulk spacetime in the SS model in a deconfined phase at finite temperature is 
given by 



ds 2 = (j^y 2 (f(u)dt 2 + StjdxW + dx 4 2 ) + {^j' 2 (u 2 dn 2 A + 



3/2 / D_ \ 3 /2 

/(«) 

The four-form field strength, the dilaton, and the curvature radius of the spacetime are 



27TiV , / U \ 3/4 ATI o 
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kD4 



respectively, where f(u) = 1 — u^/ii 3 , iiy = 167r 2 i? 3 34 T 2 /9. Xa is the compactified coordinate 
transverse to the probe D8/D8 branes with arbitrary periodicity 2ttR. The volume of the 
unit four-sphere Qa is denoted by Va and the corresponding volume 4- form by e^. Fu-) is the 
4-form field strength, l a is the string length and g s is the string coupling. 

In the SS model, the chiral symmetry dynamics is taken into account, by construction, 
in the form of the dynamics of the flavour branes, D8 and D8. The DBI action of D8-D8 is 

Sds = -/is J (fXe'^TryZ-detigMN + 2ita'F MN ) (1) 

where Fmn is the field strength of the flavour group U(Nf) on the branes. It is given by 

F = dA + %A A A. (2) 
The U(Nf) gauge field A can be decomposed into SU(Nf) part A and U(l) part A: 

A = A + —==A, (3) 

where only the diagonal U(l) will be turned on here. Lastly, gMN is the induced metric on 
the D8-branes world volume. 

In the deconfined phase, the equation of motion from the action of D8-D8 provides 3 
possible configurations: (i) connected D8-D8 without sources in the bulk representing the 
vacuum state and (ii) the parallel configuration of both D8-branes and D8 representing the 
Xs-QGP. Another stable configuration (iii) is the connected D8-D8 branes with the D4- 
brane as the baryon vertex submerged and localized in the middle of the D8 and D 80 In 



1 Actually, the quark matter, represented by the connected D8-D8 branes with radial strings stretching 
out to the horizon, is another possible configuration satisfying the equation of motion. However, it was 
found that this phase is thermodynamically unstable to density fluctuations by Bergman, Lifschytz, and 
Lippert [22]. 
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(a) (b) (c) 



Figure 1: Different configurations of D8 and D8-branes in the background field following 
the Sakai-Sugimoto model that are dual to the phases of (a) xs~QGP, (b) vacuum and (c) 
multiquark phase. 



this model, we assume that the hanging strings shrink to approximately zero and the only 
apparent strings are the k r radial strings. The three configurations are shown in Fig. [TJ We 
will consider the thermodynamic properties of only the last multiquark configuration. The 
action of the exotic multiquark phase is given by 

S = Sd8 + Sda + Sf\i (4) 

where Sds is the DBI action of the connected D8-branes, Sda represents the DBI action of 
the D4-brane wrapped on S A and Sfi is the action of k r radial strings extending from the 
baryon vertex down to the horizon. For simplicity, we ignore the distortion of the baryon 
vertex due to the Chern-Simon term j23j [23] . 

The DBI action of the D8-D8-brane coupled to the diagonal U(l) gauge field is given by 



S m = N / du U yf{u){x' A {u)f + u~ 3 (l - K(w)) 2 ), (5) 



where the constant Af = (/is^Nffl^VsR 5 ) / g s , and the rescaled U(l) diagonal field a = 
2ixa'A/(R^JW f ). Posit ion of the vertex is denoted by w c , it is determined from the equilib- 
rium condition of the D8-D4-strings configuration (see Appendix A of Ref. [17]). The source 
action of the D4 and strings, Sda + Sfi, are given by 



S source A^ d 



^m c a//(m c ) + n s (u c - ut) 



(6) 



where n s is the number of radial strings k r in the unit of N c . The number of radial strings 
n s represents the colour charges of a multiquark. For a fixed number of k^, one of the radial 
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strings can merge with another radial string from another multiquark and form a colour- 
binding potential between the two in a similar way holographic meson is formed between a 
quark and an antiquark. 

The £/b(1) symmetry corresponds to the U (l)-diagonal part of the global flavor symmetry, 
U(Nf), which is provided by the Nf flavor branes. Naturally, the baryon chemical potential, 
conjugating to the {7b (1) charge, in the gauge theory side can be identified with the boundary 
value of the zero component of the gauge field in the flavor branes, i.e. A , conjugating to 
the {7(1) "electric" charge. For convenience, our normalized baryon chemical potential is [25] 

fi = a (oo). (7) 

In gauge-gravity duality, we identify the grand canonical potential density in the gauge 
theory side in the form of the D8-branes action evaluated with the classical solution 



fi(/f) = -TjS D $[T,x' A (u),d (u)} d (8) 

With the additional source term, the free energy is in the form of the combination of the 
Legendre-transform of the grand potential and the source action, Eqn. fl6]). The baryon 
chemical potential is simply the derivative of the free energy with respect to its conjugate, 
i.e. the baryon number density, at a particular temperature: 

fJ> = g^jj [Sbs[T, x' 4 (u), d(u)]ci + S somce (d, u c )J (9) 

where the Legendre-transformed action Sds is given by 

POO 

Sr>$[T,x±(u),d(u)] = S D8 [T, x' 4 (u) , a (u)] + J\f I d(u)a' du (10) 



11 



where d(u) is the electric displacement. It is a constant of the configuration given by 
7/ \ 1 SS DB ua' (u) 

diu) = — -7 .,,. . = . = = const. (12) 

Af5a' (u) ^f{u){x' 4 {u)f + u-^l-{d' Q {u)y) 

Note that the Legendre transformation changes the dependence on the variable clq{u) in 
Sds to d(u) in Sds- As a result, the grand potential as a function of the baryon chemical 
potential is transformed into the free energy as a function of the baryon number density. 
Another constant of the configuration is 

/ / / \\2 1 r /(n)(n 8 + n 3 rf 2 ) i-i 

Mu)) = ^/>)[ J =const -' ( 13 ) 



where F is a function of u c , d, T and n s , given by 

' 9/ ( 
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where rj c = 1 + ^ + ^ n svfc- For convenience, here and henceforth, /, f c and / are 

used to represent f(u), f(u c ) and /(«o), respectively. Note that this form of F is derived 
from the force condition at the cusp u c . The detailed calculations are given in the Appendix 
of Ref. pi]. 

With the equation of motion for X4, Eqn. (1131) . and the separation between D8- and 
D8-branes L being fixed to L = 2 f°° x' A {u)du = 1, we obtain [22] 



a n (u) + 



1 



A/" dd 



(15) 



T,L ,u c 



where the second term is the contribution from the sources, fi SOU rce- From these relations, we 
can then study the thermodynamic properties of the multiquark phase. The phase diagram 
of the multiquark nuclear phase is studied in Ref. p2] when the colour-binding interaction 
is neglected. It is found that multiquarks are preferred thermodynamically over the other 
gluon-deconfined phases for the large density and intermediate temperature below the chiral 
symmetry restoration temperature. 



3 Calculations of the equation of state 



Thermodynamic properties of the nuclear /exotic matter phase can be described by the equa- 
tion of state. First, we will investigate the relations between the pressure and the number 
density. From the previous section (see also Ref. [H]), the grand potential density and the 
chemical potential of the nuclear/exotic matters are given by 



n 



du 



du 



-1/2 



u ' 



f(u)(u 8 + u s d 2 )_ 
F 2 



Vu 5 + d 2 ' 
- 1/2 d 



f(u)(u 8 + u 3 d 2 )\ V^Td^ 3 



o u cV f( u c) + n s (u 



(16) 

-u T ) (17) 



respectively. 

Since the differential of the grand potential Gq can be written as 



dG< 



-PdV - SdT - Ndfi 



(18) 



where the state parameters describing the system P, V, S, T, N are the pressure, volume, 
entropy, temperature, the total number of particles of the system respectively. Since the 
change of volume is not our concern, we define the volume density of Gq, S and N to be Q, 
s and d, respectively. Therefore, we have, at a particular T and /x, 



P=-G a /V=-n(T, t j). 
By assuming that the multi-quark states are spatially uniform, we obtain 

dP 



(19) 



(20) 
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Using the chain rule, 



dP 
~dd 



T 



dd 



d. 



so that 



f d 

P(d,T,n s ) = fi(d,T,n s ) d — I f/,(d',T,n s ) d(d') 

Jo 



(21) 



(22) 



where we have assumed that the regulated pressure is zero when there is no nuclear matter, 
i.e. d = 0. 

In the limit of very small d, u c approaches Uq, Vc becomes t]q + 0(d), where rjo is defined 
to be r] c with u c replaced by uq. From Eqn. ffTTj) . the baryon chemical potential can then be 
approximated to be 

<h i°< (i - £ - 3) " /2 



1^ /^source 



d< 



du 



«- 5/2 (l-^) I (23) 



where fi source = \ u c\/ f(u c ) + n s (u c — u?), and we have neglected the higher order terms of 
d. By using the binomial expansion, the above equation becomes 



V - Use 



d< 



u 



-5/2 



du 



/« 8 



foul 



1 



fu* - f ul V 9/o « 5 



1 \ d 2 



w 



a d — f3 (n s )d 3 



where 



a 



/3o(n s 



du- 
1 

o 

du 



u 



-5/2 



U 



fov% 
5/2 



fou 3 



o./l f° u o \f u8 ~ fo U V 9 /o U l 



u: 



(24) 

(25) 
(26) 



By substituting Eqn. (I24p into Eqn.f l22p . we can determine the pressure in the limit of very 
small d as 



P ~ - 3 ^ Us h 4 
2 4 

In the limit of very large d and relatively small T, 



/* /^source 



du 



du 



f c u 6 c u c + a 



9/c 



fu* 



U 5 + d 2 



d 



d 



Vu 5 + d 2 2 



+ -u*J c d< 1 - 



ft_ 

9/c 



Vn 5 + d 2 



du 



d 



fu 3 {u 5 + d 2 fl 2 



^ 2/5 r(|)r(A) ^ / _ ^ /5 r(-§)r(j§) 
5 r(|) + io I 9fJ d r(l) 



(I) 



(27) 

(28) 
(29) 
(30) 



where we have used the fact that the lower limit of integration u^/d 2 is approximately zero 
as d is very large. Again by using Eqn. (j22p . we obtain 

p -l( £ f5P) < ' 7/5 (3i) 

And the energy density can then be found via the relation dp = fid(d). 

Next we consider the entropy of the multiquarks phase. From the differential of the free 
energy, 

dF E = -PdV - SdT + pdN, (32) 

the entropy is given by 

S = (33) 



The entropy density can then be written as 



dTE (34) 



dT 



where Te is the free energy density which relates to the grand potential density as Te 
Q + fid. Since we have the pressure P = —Q, we can write 



dP_ _ fd\i_ 
dT ~ \&T 



d. (35) 



For both small d and large d, we can see from the formula of the pressure (see Eqn. (I27D . (I3~I1) . 
noting that a , (3 is insensitive to temperature) and the chemical potential (see Eqn. (1241) . (130|) ). 
that the dominant contribution comes only from p SOU rce, thus 



djjL s 



dT 



d. (36) 



The baryon chemical potential from the D8-branes is insensitive to the changes of temper- 
ature. This implies that the main contribution to the entropy density of the multiquark 
nuclear phase comes from the source term namely the vertex and strings. 
Since 



-u c yj f(u c ) + n s (u c - u T ) ) , (37) 



dpsource 9 

dT df \3 

d/Wce (^) T ' 



dT r / \3 9 



u 2 Q 



n s — — , (38) 



where we have used the fact that u c is approximately constant with respect to the tempera- 
ture in the range between the gluon deconfinement and the chiral symmetry restoration (see 
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Figure 2: The graphs show the relations between u c and T at small density (left) and at 
large density (right). 

Fig. [2]). Therefore, we obtain 



For small n s , the entropy density is proportional to T 5 . When n s gets larger (carrying 
colour charge), the entropy density becomes dominated by the colour term s oc n s T. This 
is confirmed numerically in Section 4. It has been found that the entropy density of the 
Xs-QGP scales as T 6 [22] corresponding to the fluid of mostly free quarks and gluons. We 
can see that the effect of the colour charge of the multiquarks as quasi-particles is to make 
them less like free particles with the temperature dependence ~ n s T, i.e. much less sensitive 
to the temperature. 

It is interesting to compare the dependence of pressure on the number density, Eqn. ( 1271) 
and ( 13 lj) . to the confined case at zero temperature studied in Ref. [28]. The power-law rela- 
tions for both small and large density of the confined and deconfined multiquark phases are 
in the same form (for n s = 0). The reason is that the main contributions to the pressure 
for both phases are given by the D8-branes parts and they have similar dependence on the 
density for both phases. For the deconfined multiquark phase, the additional contributions 
from the source terms in Eqn. ( fl71) . fi SOU rce, are mostly constant with respect to the den- 
sity (this is because u c becomes approximately independent of d for small and large d limits). 




(39) 
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Consequently, when we substitute into Eqn. (j22p . the constant contributions cancel out and 
affect nothing on the pressure. 

On the contrary, the entropy density for the deconfined phase is dominated by the contri- 
butions from the sources namely the vertex and strings. The contribution of the D8-branes 
is insensitive to the change of temperature and therefore does not affect the entropy density 
significantly. The additional source terms, however, depend on the temperature and thus 
contribute dominantly to the entropy density. Once the temperature rises beyond the gluon- 
deconfined temperature, entropy density will rise abruptly (for sufficiently large density d) 
and become sensitive to the temperature according to Eqn. (139|) . due to the release of quarks 
from colourless confinement appearing as the sources. However, we will see later on using 
the numerical study in Section 4 that for low densities and for small n s , the numerical value 
of the entropy density is yet relatively small. 

4 Numerical studies of the thermodynamic relations 

From the analytic approximations in the previous section, we expect the pressure to appear 
as straight line in the logarithmic scale for small and large d with the slope approximately 2 
and 7/5 respectively. The relation between pressure and density of the multiquarks from the 
full expressions can be plotted numerically as are shown in Fig. [3]l5j The pressure does not 
really depend much on the temperature and we therefore present only the plots at T = 0.03. 
Remarkably, the transition from small to large d is clearly visible in the logarithmic-scale 
plots. The transition occurs around d c ~ 0.072. Interestingly, as is shown in Fig. the 
multiquarks with larger n s has lower pressure than the ones with smaller n s for d < d c 
and vice versa. The dependence on n s remains to be seen for small d as we can see from 
Eqn. ( 127]) . For large d, the n s -dependence is highly suppressed as predicted by Eqn. (l3Tj) . 

The entropy density as a function of the temperature for various ranges of the density is 
shown in Fig. El The temperature dependence for both small and large d are the same, ~ T 5 
at the leading order. The ci-dependence is linear and thus appears as separation of straight 
lines in the logarithmic-scale plot. For n s > 0, we can see from Eqn. (139]) that the linear 
term in T should become increasingly important. This is confirmed numerically as is shown 
in Fig. [6j The slope of the graph between the entropy density s and T in the double-log scale 
for n s = (the left plot) and n s = 0.3 (the right plot) is approximately 5 and 1 respectively. 
Regardless of the temperature dependence, it should be noted that the numerical value of 
the entropy density for small densities and low n s in Fig. [6] is quite small. 

Lastly, the relations between the baryon number density and chemical potential are shown 
in Fig. [7J Temperature has very small effect on these curves and negligible for the range 
of temperature between the gluon deconfinement and the chiral-symmetry restoration. The 
baryon chemical potential depends linearly on the number density for small d. For large d, the 
relation between the chemical potential and number density becomes [i ~ d 2 ^ '. Interestingly, 
the multiquark quasi-particles behave more like fermions as a result of being the electric 
response of the DBI action [22] . 
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Figure 3: Pressure and density in logarithmic scale at T = 0.03. 
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Figure 6: Entropy and temperature in logarithmic scale for n s = (left), 0.3 (right). 
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Figure 7: The baryon chemical potential and number density in logarithmic scale at T = 0.03. 

5 Gravitational stability of the dense multiquark star 

When a dying star collapses under its own gravity, it is generically believed that the degen- 
eracy pressure of either electrons or neutrons would be able to stop the collapse to form a 
white dwarf or a neutron star. If the star is more massive than the upper mass limit of the 
neutron star, it would collapse into a black hole eventually. The mass limit of the neutron 
star is sensitive to the physics of warm dense nuclear matter but little is known about the 
equation of state of nuclear matter under high temperature and large density. Even though 
the original mass limit of the neutron star estimated by Oppenheimer, and Volkoff was only 
0.7 solar mass [2T], the new limit when the nuclear interactions are included could be as large 
as 2.5 solar mass [27]. Under extreme pressure and density, the quarks within hadrons could 
be freed and wander around the interior of the star. In other words, quarks are effectively 
deconfined from the localized hadrons but confined by gravity within the star. Using the bag 
model to describe the state of being confined by gravity but possibly deconfined from the 
hadrons, it turns out that quark matter phase, e.g. strange star, is possible under extreme 
pressure and density. 

However, physics of the deconfinement is largely unknown due to the non-perturbative 
nature of the strong interaction and the difficulty of lattice approach to deal with finite 
baryon density situation. The bag model are not always served as a reliable theoretical tool 
to explore the behaviour of quarks in the dense star when the deconfinement exists. It is 
therefore interesting to use the equation of state of the deconfined nuclear matter from the 
holographic model to investigate the behaviour of the dense star as a complementary tool 
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to the bag model and other approaches. 

In this section, we will consider a hypothetical multiquark star containing only the mul- 
tiquark matter with uniform constant temperature. The relations between pressure and 
density will be adopted directly from the holographic model as the equations of state of the 
quasi-particles. Since the pressure and density have very small temperature dependence for 
the range of temperatures under consideration, the results are valid generically. 

A study into the gravitational stability of a spherically symmetric dense star can be 
performed using the Tolman-Oppenheimer-Volkoff equation [21]. It is known that the spher- 
ically symmetric dense star has metric in a form 

ds 2 = A{r)dt 2 - B{r)dr 2 - r 2 dQ 2 . (40) 

After substituting into the Einstein field equation, we obtain the following relations, 

B(r) = (l-^P-Y, (41) 



dA*(r) 



dr 



8npr 2 , (42) 



and 

dP{r) (p + P)A'(r) (p + P) 8?rPr 3 + A*(r) 
dr ~ 2 A(r) ~ 2 r(r - A*(r)) ^ ' 

The last equation is known as the Tolman-Oppenheimer-Volkoff (TOV) equation. The accu- 
mulated mass of the star, M(r), is given by A*(r) = 2M(r). It has been shown in Ref. [29] 
that the chemical potential can be defined through the background metric in the form of 

p(r) = pH — It will automatically solve the TOV equation. Note that the constant ep is 
V A ( r ) 



arbitrary. Since 
the TOV equation becomes 



Ar)_ l A'jr) 

p(r) 2 A(r) ' 1 ' 

d -^= iP+P m (45) 

dr p[r) 

Together with the first law of thermodynamics p + P = pd, the TOV equation then takes 
the following form, 

1 fdP\ 

^ = 1 d ( d )- ( 46 ) 



d \dd _ 

Obviously, the chemical potential can be determined, as a function of the number density: 

f d 1 fdP\ 

M°0 = / -\-Zr) dr l + Vonset, (47) 

Jo V V 9rj J 

where p onS et = ^(d = 0). Additionally, considering from the TOV equation together with 
the first law of thermodynamics, the density dp = pd(d) can be integrated to 



dr). (48) 
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For a power-law equation of state, P = kd x , the chemical potential, Eqn. (I47p . becomes 

p(d) = yTJ^ 1 + /i °" set ' ( 49 ) 
and eventually the equation of state is given by 

1 fP\ 1/x 

P = T T^P + Ponset [ ~T ] ■ (50) 



— 1 \ K 

In our holographic model of multiquarks, the relation between pressure and density has a 
unique power-law behaviour, as is also found in Ref. [22] for the case of normal baryon (n s = 
0). This is shown in Fig. [3B For small d, P oc d 2 (n s = 0) and for large d, P oc d 7 ^ 5 . 
The dependence on n s becomes significant when the density d is small and the equation of 
state can be approximated by P ~ ad 2 + (3d 4 . Since there are two power-laws governing, we 
need to match the solutions from the two regions together (i.e. core and crust). The number 
density where the equation of state changes from the large-ci to the small-d is denoted by d c . 

For n s = 0, at the transition point d = d c , the energy density is given by Eqn. (1501) . 



Pc — -TZ r + Ponsetdc, (51) 

where P = k'd x ' (Eqn. ( 1271) suggests that A' = 2) is the equation of state of the small d 
region. We recalculate the relation Eqn. (|47j) . (148]) for the large d region which match with 
this p c to be 

p = ^ + \ k (-—--J—\ (52) 



p = pc + yz\ p + Pc 



k) 



d. 



Numerical results and Eqn. ( l31~j) suggest that A = 7/5 for the large d region. 

For n s > 0, assume the equation of state for small d is in the form of P = ad Xl + 
bd X2 (Eqn. ( 1271) suggests that Ai^ = 2,4), the chemical potential and energy density for the 
small d region become 

Axac^- 1 \ 2 bd X2 ~ 1 , , 

/i = Ponset + — ; — + -: — , (54) 

Al — 1 A2 — 1 



ad Xl bd X2 

P = Ponsetd+- - + t -. (55) 

A{ — 1 A2 — 1 



We obtain the transition density in the similar fashion, 



ddc 1 bd c 2 rr^ 
Pc = Ponsetdc + T T + T (56) 

Ai — 1 A2 — 1 
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Numerical results show that for large d, the effect of n s is negligible. Therefore, the 
baryon chemical potential and the density for the large d region are again given by Eqn. (j52p 
and f )53|) . The equations of state, Eqn. f l50|) .f l53|) as well as the corresponding relations for 
n s > case, are in the mixed form containing both the quasi-particle nonlinear terms and 
the linear term. The linear term is roughly pii near ~ 2.5P and the quasi-particle term is 
approximately p quasi « P 5/7 . 

We can solve the TOV equation when the equations of state are given as above by starting 
from the core of the star out to the surface. As we go from the center towards the surface of 
the star, the density decreases until it reaches a critical value p c . This density corresponds 
to the number density d c where the power-law changes from P ~ d 7 ^ to P ~ d 2 (see Fig. [Si- 
ll]). For the crust region where the density p < p c , multiquarks obey a different equation 
of state given by Eqn. (1501) . The radius of the core is defined to be the distance R core 
where p(R CO re) = pc and the surface of the star is defined to be the radial distance R where 
p(R) = 0. 
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Figure 8: The relation between mass and central density of the multiquark star for multi- 
quarks with n s = (upper), 0.3 (lower). 

For n s = 0, numerical fittings suggest that k = 10 -0 4 , A = 7/5, d c = 0.215443, p c = 
0.564374 (core) and k! = 1, A' = 2, p onset = 0.17495 (crust). For n s = 0.3, good fit parameters 
are k = 10" a4 , A = 7/5, d c = 0.086666, p c = 0.490069 (core) and a,b = 0.375, 180.0, ; X 1>2 = 
2,4; p onse t = 0.32767 (crust). Varying the central density p of the star, we obtain the mass- 
density relation in Fig. [HJ Each curve has two maxima, a larger one in the small density 
region and a smaller one in the large density region. Each maximum corresponds to each 
power-law of the equation of state, the low density to the crust and the large density to 
the core. Interestingly, the contribution to the total mass of the multiquark star comes 
dominantly from the crust. This is shown in Fig. |9j Even though the density is much lower, 
the volume of the crust is proportional to the second power of the radius and thus makes 
the contribution of the crust to the total mass larger than the core's. 
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Figure 9: The accumulated mass distribution in the hypothetical multiquark star for the 
central density po = 20 and n s = 0. The inner (outer) red (dashed-blue) line represents the 
core (crust) region. 





Figure 10: The density, and pressure distribution in the hypothetical multiquark star for the 
central density p — 20 and n s = 0. The inner (outer) red (dashed-blue) line represents the 
core (crust) region. 
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Figure 11: The adiabatic index at constant entropy (r) and the sound speed (c s ) distribution 
in the hypothetical multiquark star for the central density p = 20 and n s = 0. The 
inner (outer) red (dashed-blue) line represents the core (crust) region. 
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Figure 12: Comparison of the accumulated mass distribution in the hypothetical multiquark 
star for the central density po = 20 between n s = and 0.3. The (dashed) blue line represents 
the crust region of multiquark star with n s = 0.3 (0). The red lines represent the core region 
of which both cases are almost the same. 
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Figure 13: Comparison of the density, and pressure distribution in the hypothetical multi- 
quark star for the central density po — 20 between n s = and 0.3. The (dashed) blue line 
represents the crust region of multiquark star with n s = 0.3 (0). The red lines represent the 
core region of which both cases are almost the same. 
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Figure 14: Comparison of the baryon chemical distributions in the hypothetical multiquark 
star for the central density p = 20 between n s = (left) and 0.3 (right). The solid (dashed) 
red (blue) line represents the core (crust) region. 
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Figure 15: The adiabatic index at constant entropy (r) and the sound speed (c s ) distribution 
in the hypothetical multiquark star for the central density po = 20 and n s = 0.3. The 
inner (outer) red (dashed-blue) line represents the core (crust) region. 
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Figure [TU] shows the pressure and density distribution within the multiquark star for the 
case of n s = for the central density p = 20. Even though the density and pressure decrease 
rapidly with respect to the radius of the star, they never quite reach zero. It turns out that 
when the density and pressure reach the critical values where the equation of state changes 
into the different power-law for small d, the crust region continues for a large fraction of the 
total radius of the star. This makes the crust mass contribution to the total mass of the star 
dominant as is shown in Fig. O 

Some remarks should be made regarding the hydrodynamic properties of the multiquark 
phase (taken as nuclear liquid). At constant temperature and entropy, we can define the 
adiabatic index 



r = 



p dP 
P~dp'' 



^c 2 
P s 



(57) 
(58) 



where c s is the sound speed in the multiquark liquid. They depend on the equation of state 
of the multiquark and their distributions within the multiquark star are shown in Fig. [11] 
for n s = 0. The sound speed never exceeds the speed of light in vacuum. It is also found 
that the adiabatic index and the sound speed change within a small fraction as the central 
densities are varied for a given n s . 

The multiquark star with n s = 0.3 (having colour charges) converge to a smaller mass 
and radius at high central density (Fig. [T2|) . Multiquarks with colour charges has lower 
pressure (and therefore smaller density) than the colourless ones for small density (Fig. [T3|) . 
This smaller pressure makes the coloured multiquark star smaller and thus less massive than 
the colourless one. In more realistic situations, all of the possible multiquarks with varying 
n s coexist in the multiquark phase. The mass limit and mass radius relation will vary 
between the two typical cases we consider here. Since the equations of state are found NOT 
to be sensitive to the temperature within the range between the gluon deconfinement and 
the chiral symmetry restoration, our results should also be valid even when the temperature 
varies within the star (but not too high and too low, of course). 
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Figure 16: The relation between mass and radius of the multiquark star with (a) n s = 0, (b) 
n s = (red) and n s = 0.3 (black). 



21 




Figure 17: The relation between mass and radius of the core of the multiquark star with (a) 
n s = 0, (b) n s = (blue) and n s = 0.3 (black). 

The baryon chemical potential distributions in the multiquark star for n s = 0, 0.3 are 
shown in Fig. [TH In the core region, the chemical potential distributions of both cases are 
similar due to the similarity of the equations of state for large density. A small jump of the 
chemical potential at the transition radius between core and crust region is the artifact from 
the power-law approximation. The value of the chemical potential at the transition radius 
from the full expression which we used in the numerical simulations is slightly different from 
the approximated value using the power-law. 

The adiabatic index and sound speed of the multiquark phase for n s = 0.3 are shown in 
Fig. [15] The adiabatic index is higher than n s = case but the sound speed in the low density 
region is distinctively higher. Around the transition density, the sound speed reaches the 
maximum value of about 0.986 of the speed of light in vacuum. For both n s = 0, 0.3 cases, it 
is obvious that the adiabatic index is closer to 1 in the core reflecting the fact that the density 
distribution is more condensed in the core region. The adiabatic index reaches A' = 2 at the 
star surface since the the equation of state at zero density is P oc p x (i.e. T(p — > 0) = A' 
for Eqn. (|50l). 

The spiral relation between mass and radius of the multiquark star is shown in Fig. [16j 
As the central density is increasing, the mass and radius of the n s = (0.3) multiquark 
star converge to the value of 0.659 (0.440) and 3.132 (1.704) respectively. For the core, the 
mass and radius of the core for n s = (0.3) converge to the value of 0.108 (0.169) and 
0.471 (0.737). 

Finally, we would like to estimate these limits of mass and radius in the physical units. 
Since our dimensionless quantities are related to the physical quantities through conversion 
factors given in Table [TJ (Appendix A), both physical mass and radius vary with the energy 
density of the nuclear matter phase as oc 1/ -^/energy density scale. For a multiquark nuclear 
phase with energy density scale 10 GeV/fm 3 , the conversion factor of the mass and radius 
are 5.91M<, ; ar and 8.71 km respectively. This would correspond to the converging mass and 
radius (in the limit of very large central density) of 3.89 (2.60)M so i ar and 27.29 (14.85) km 
for n s = (0.3) multiquark star respectively. 

In realistic situation, the nuclear phase in the outer region could lose heat out to the 
space in the form of radiation. The nuclear matter in the outer region of the crust will cool 
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down and mostly become confined into neutrons and hadrons (e.g. hyperons, pions). This 
would make the multiquark crust to end at shorter radius than the estimated value and 
render the multiquark star to be smaller and less massive than the estimated values in the 
hypothetical prototype. For example, for the energy density scale 10 GeV/fm 3 , the critical 
density is p c ~ 1.5 x 10 18 kg/m 3 (n s = 0). This is still a sufficiently large density for the 
neutron layer to be formed. If the temperature of the nuclear matter in the crust region 
falls below the deconfinement temperature, the multiquarks will be confined into extremely 
dense neutrons and hadrons instead. For a typical neutron star, the distance of the neutron 
layer out to the star surface is roughly 5-6 km [30]. If we add this number to the radius of 
the multiquark core, 0.471 x 8.71 ~ 4.10 km, we end up with a more realistic estimation 
for the multiquark star with radius ~ 10 km. Regardless of the name, only the core region 
is in the deconfined multiquark phase and the content of the outer layers are the confined 
nucleons. 

6 Discussions and conclusion 

In the gluon-deconfined phase of the general Sakai-Sugimoto model, multiquark states can 
exist in the intermediate temperatures below the chiral symmetry restoration temperature 
provided that the density is sufficiently large. They are stable and preferred thermodynam- 
ically over other phases and thus they can play an important role in the physics of compact 
warm stars. By analytic and numerical methods, we demonstrate that the equation of state 
of the multiquark nuclear matter can be approximated by two power-laws in the small and 
large density region. Roughly speaking, the pressure is proportional to d 2 and d 7 ^ 5 for the 
small and large number density (d) regions respectively. 

It is also found that the effect of the colour charges of the multiquark is to reduce the 
pressure of the multiquarks when the density is small. At higher densities, multiquarks 
with colour charges exert slightly larger pressure than the colourless ones. The temperature 
dependence of the entropy density shows an s oc T 5 relation and the colour charge dependence 
s colour oc n s T (see Fig. [6] and Eqn. fl39|) ). This implies that the multiquarks with colour 
charges have larger entropy but their number of degrees of freedom depend less sensitively 
on the temperature. Multiquarks in the deconfined phase behave like quasi-particles with 
the entropy density being less sensitive to the temperature than the gas of mostly free gluons 
and quarks in the Xs-QGP phase. 

Using the power-law equations of state for both small and large density regions, a spheri- 
cally symmetric Einstein field equation is solved to obtain the Tolman-Oppenheimer-Volkoff 
equation. By solving this equation numerically, we establish the mass, density and pressure 
distribution of the hypothetical multiquark star. It turns out that the multiquark star is 
separated into two layers, a core with higher density and a crust with lower density. Mass 
limit curve is also obtained as well as the mass sequence plot between the mass and radius 
of the multiquark star. They show typical spiral behaviour of the star sequence plots. The 
mass limit curve shows two peaks corresponding to the equation of state of the small and 
the large density. Analyses show that the most contribution to the total mass is mainly 
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from the crust. The adiabatic index at constant entropy, T, and the sound speed, c s , of the 
multiquark nuclear phase within the star are calculated numerically. For large density, T is 
approximately close to 1 and c s is roughly within range 0.6 — 0.7 of the speed of light. For 
small density, V is in the range 1.3 — 2.0 (2.0 — 3.0) and c s is roughly — 0.85 (0 — 0.99) for 
multiquark with n s = (0.3). 
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A Dimensional translation table 



quantity 



dimensionless variable 



physical variable 



pressure 



density 



mass 



P 

M 



P 




radius 



r 



r r 



Table 1: Dimensional translation table of relevant physical quantities, r$ 
(^■(energy density scale)) 
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